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摘 要 : 科学 认识 生境 质量 的 时 空 变化 特征 及 其 发 展 趋势 ,对 黄河 中 游 生态 保护 修复 以 及 高 质量 发 展 具有 重要 的 
现实 意义 。 本 文 基于 1992 一 2020 年 欧洲 航空 航天 局 逐年 的 土地 覆盖 数据 ,采用 InVEST 模 型 估算 生境 质量 ,运用 
Sen+MK 检验 和 Hurst 指 数 分 析 方 法 ,探究 了 退耕 还 林 ( 草 ) 工 程 实施 前 后 黄河 中 游 的 生境 质量 时 空 变 化 和 发 展 趋 


势 。 研 究 结果 表明 :(1) 整体 上 看 ,29 a 来 黄河 中 游 地 区 生境 质量 整体 较 好 ,平均 生境 质量 指数 为 0.652 ,但 呈现 小 
幅 下 降 趋 势 (每 10a 下 降 0.004), 从 局 部 上 看 ,退耕 还 林 ( 草 ) 区 生境 质量 呈现 上 升 趋势 (每 10 a 上 升 0.174) ,其 中 ， 
2011 一 2020 年 上 升 趋势 最 为 显著 ,10 a 增加 了 0.656;(2) 2001 一 2020 年 黄河 中 游 地 区 生境 质量 发 生 时 空 变 化 的 趋 
势 类 型 以 明显 降低 为 主 ,但 从 不 同时 间 段 的 发 展 来 看 ,2011 一 2020 年 较 2001 一 2010 年 生境 质量 向 好 发 展 的 面积 增 
多 ,体现 在 退耕 还 林 ( 草 ) 区 的 生境 质量 上 升 ;(3) 未 来 黄河 中 游 地 区 生境 质量 发 生 时 空 变化 趋势 类 型 以 明显 降低 和 


明显 上 升 为 主 , 且 可 持续 性 较 强 ,明显 降低 区 主要 分 布 在 汾 渭 平原 以 及 城市 集中 发 展区 域 ,明显 上 升 区 主要 位 于 山 


西 省 .陕西 省 以 及 内 蒙古 自治 区 三 地 交界 片区 。 
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黄河 流域 是 中 国 重要 的 生态 屏障 ,对 维护 区 域 
生态 安全 具有 十 分 重要 的 作用 "。2019 年 9 月 18 
日 ,习近平 总 书记 在 黄河 流域 生态 保护 和 高 质量 发 
展 座谈 会 上 提出 “黄河 流域 生态 保护 和 高 质量 发 展 
上 升 为 国家 重大 战略 ” ,在 大 力 推进 生态 保护 修复 
的 黄河 流域 中 游 显得 尤为 重要 ” , 随 着 黄河 中 游人 
类 改造 自然 界 的 规模 和 力度 不 断 增强 ,土地 利用 强 
度 及 方式 的 时 空 变化 使 得 生态 系统 结构 与 格局 改 
变 , 对 生境 要 素 .生境 区 域 之 间 的 物质 能 量 循环 产 

影响 汪 , 即 改 变 了 生境 质量 。 退 耕 还 林 ( 草 ) 工 程 已 
实施 二 十 几 年 ,黄河 中 游 土 地 利用 转变 格局 明显 ,从 
而 改变 了 区 域 生境 质量 。 因 此 ,科学 认识 生境 质 
量 的 时 空 变化 特征 及 其 发 展 趋势 ,对 黄河 中 游 生态 
保护 修复 治理 以 及 高 质量 发 展 具 有 重要 的 现实 


近年 来 ,学 者 们 广泛 应 用 生态 系统 服务 评估 与 
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权衡 模型 (Integrated Valuation of Ecosystem Services 
and Trade-offs,InVEST) 中 的 生境 质量 模块 (Habitat 
Quality,HQ) 来 估算 生境 质量 ,并 具有 较 好 的 效果 "。 
土地 利用 / 获 盖 数据 是 生境 质量 测算 的 核心 数据 ， 
中 科 院 资源 环境 科学 数据 中 心 可 提供 5a 时 间 尺 度 
的 数据 (分 状 率 30 m)" ,国家 基础 地 理 信息 中 心 
能 够 提供 10 a 时 间 尺 度 的 地 表 履 被 产品 (分 辩 率 
30 m) ,学 者 们 也 可 根据 研究 需要 ,下 载 处 理 多 
时 期 遥感 影像 ,采用 目 视 解 译 或 面向 对 象 分 类 等 方 
法 解 译 出 各 时 期 的 土地 利用 类 型 “" ,在 现 有 土地 利 
用 类 型 数据 基础 上 ,运用 CA-Markov 模 型 模拟 预测 
未 来 时 期 的 土地 利用 类 型 "中 。 在 年 际 时 间 尺 度 
上 ,可 采用 NASA 提供 的 MODIS MCD12Q1( 分 辩 率 
500m) 土 地 利用 /覆盖 产品 数据 ”。 然 而 ,中 高 空间 
分 辨 率 土 地 利用 /覆盖 数据 时 间 间 隔 较 长 ,时 间 序 列 
数据 难以 获取 ,满足 时 序 条 件 的 土地 利用 / 窗 盖 数据 
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空间 分 辨 率 相对 较 低 ,数据 应 用 上 限制 了 长 时 间 序 
列 生 境 质量 时 空 变化 规律 分 析 。 而 由 欧洲 空间 局 
(European Space Agency, ESA) ii] 备 的 (Climate 
Change Initiative-land Cover, CCI-LC ) Hu R #2 iti 28 AY 
产品 ,分 类 精细 时间 序列 最 长 (1992 一 2020 年 )、 空 
间 分 辨 率 较 高 (分 辨 率 300 m) ,在 探究 退耕 还 林 
( 草 ) 工 程 实施 前 后 黄河 中 游 生境 质量 时 空 变化 特 
征 具 有 明显 的 优势 。 

学 者 们 通过 生境 质量 时 空 变化 分 析 , 在 土地 整 
治 和 生态 保护 恢复 等 监测 评价 方面 具有 较 好 的 应 
用 。 如 在 土地 整治 前 、 中 、 后 3 个 时 期 ,生境 质量 经 
历 了 先 下 降 后 上 升 的 过 程 ”, 农 牧 交错 带 在 开展 一 
系列 生态 保护 和 生态 恢复 工程 后 ,区 域 生境 质量 总 
体 趋 于 好 转 汪 ,退耕 还 林 ( 章 ) 政 策 实 施 可 促使 生境 
质量 转 好 ”3。 同 时 ,不 合理 的 人 类 活动 会 导致 生境 
质量 退化 ,如 在 黄河 湿地 自然 保护 区 ,建设 用 地 面 
积 的 增加 和 湿地 的 波动 减少 是 导致 生境 质量 退化 
的 主要 原因 ,黄河 三 角 洲 高 效 生态 经 济 区 的 土地 
利用 程度 与 生境 质量 呈 显 著 负 相关 性 , 且 相 关 性 不 
断 增 强 ,人 类 活动 与 生境 质量 之 间 存 在 显著 的 空 
间 依 赖 性 ,二 者 呈现 显著 的 空间 负 相 关 光 。 黄 河中 
游 地 区 生态 环境 脆弱 ,经 济 与 生态 矛盾 突出 ,在 
重点 实施 退耕 还 林 ( 草 ) 生 态 恢 复工 程 背 景 下 ,29 a 
年 来 黄河 中 游 地 区 及 其 退耕 还 林 ( 草 ) 区 的 生境 质 
量 年 际 变化 特征 和 时 空 变 化 趋势 如 何 ” 哪 些 地 区 
的 生境 质量 具备 上 升 潜力 ,哪些 地 区 未 来 需要 重点 
关注 ? 这 些 问 题 尚 不 明确 ,需要 结合 遥感 长 时 序数 
据 与 地 理 信息 空间 分 析 技 术 进 行 定量 研究 与 评价 。 


图 例 

黄河 流域 
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鉴于 此 ,本 研究 在 退耕 还 林 ( 草 ) 工 程 实施 前 
后 ,基于 1992 一 2020 年 未 年 土地 覆盖 数据 ,借助 In- 
VEST 模 型 估算 生境 质量 ,运用 Theil-Sen Median 趋 
势 分 析 方 法 ,揭示 黄河 中 游 生境 质量 时 空 演变 特 
征 , 并 采用 Hurst 指 数 探 究 未 来 生境 质量 的 发 展 持 
续 性 。 研 究 结果 可 为 黄河 流域 生态 保护 和 高 质量 
发 展 提供 科学 依据 。 


1 数据 和 方法 


1.1 研究 区 概况 

黄河 流域 中 游 地 区 是 指 内 蒙古 自治 区 呼 和 浩 
特 市 的 河口 镇 至 河南 省 郑州 市 的 桃花 峪 区 间 呈 ,地 
BS 103°53'~113°41'E,33°39' ~40°35' N, W RAK 
1234 km, 流 域 面 积 为 3.44x10 km’, FIR TE 0~3929 
m, SAS CRS AE, PEFR YA BRA BY Be E E Bae A a E 
交错 的 黄土 高 原 ( 图 1)。 该 区 处 于 我 国 半 湿 润 区 向 
干旱 区 过 渡 地 带 , 降 水 具有 很 强 的 空间 变异 性 , 空 
间 格 局 由 东南 向 西北 .由 山地 向 平地 递减 ,降水 集 
中 在 6 一 9 月 ,多 以 暴雨 为 主 , 水 土 流失 严重 。 为 了 
改善 生态 环境 ,黄河 中 游 退 耕 还 林 ( 草 ) 项 目 大 规模 
实施 ,因地制宜 的 进行 植树 造林 .恢复 植被 ,生态 环 
境 得 到 较 好 的 恢复 ,土地 利用 / 履 盖 状 况 发 生 了 巨大 
变化 ”。 
1.2 数据 来 源 与 处 理 

文中 使 用 的 土地 覆盖 数据 来 源 于 欧洲 航空 航 
天 局 研发 的 逐年 土地 覆盖 数据 集 , 其 是 目前 时 间 序 
列 最 长 (1992 一 2020 年 )、 空 间 分 辩 率 较 高 (300 mx 
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300 m) 的 全 球 土地 覆盖 产品 ,非常 适合 基于 土地 覆 
盖 变 化 的 长 时 间 序 列 、 大 范围 的 时 空 变化 研究 ”。 
从 欧 空 局 网 站 (http://maps.elie.ucl.ac.be/CCl/viewer/) 
下 载 全 球 的 NetCDF4 格 式 数据 ,人 然后 利用 QGIS 软件 
处 理 成 证 格式 数据 ,坐标 系统 为 双 CS1984, 最 后 通 
过 裁剪 处 理 得 到 研究 区 1992—2020 F RAE LH 
盖 数 据 。 该 数据 集 总 体 精 度 为 71.98% ,Kappa 系 数 
为 0.650250。 该 产品 依照 联合 国 粮农 组 织 开 发 的 土 
地 履 盖 分 类 系统 标准 分 为 22 个 大 类 ,36 个 精细 分 
类 。 为 了 便于 黄河 中 游 地 区 生境 质量 的 评估 ,基于 
精细 分 类 的 36 个 类 型 ,依据 CCLLC 土地 利用 /覆盖 
类 型 重 分 类 规则 归并 为 耕地 、 林 地 .草地 湿地 水 
体 、 人 类 造 地 表 和 未 利用 地 7 类 加 ,得 到 黄河 中 游 地 
区 的 七 大 类 土地 覆盖 类 型 。 在 本 研究 中 所 涉及 的 
退耕 还 林 ( 草 ) 区 是 2000 一 2020 年 由 耕地 转化 为 林 
地 或 草地 的 栅 格 空间 分 布 区 。 

1.3 研究 方法 

1.3.1 InVEST 模 型 ”利用 1992 一 2020 年 土地 利用 / 
履 盖 类 型 机 格 数据 ,依托 InVEST 模 型 中 的 Habitat 
Quality 模块 计算 每 年 的 生境 质量 指数 了 。 首 先 , 确 
定 生境 威胁 因子 敏感 度 和 外 界 威胁 强度 ,计算 出 生 
境 质 量 的 退化 程度 ,然后 设置 半 饱 和 常数 (4) 得 到 生 
境 质量 。 具 体 公 式 如 下 : 


d 
-| | (线性 衰退 ) 
i = rmax (1 ) 
ef {32°}, (指数 衰退 ) 
>>> er xr Xi XB, XS, (2) 
J 7 7 Se, ki ‘) J 


De 
on o 


SUP tiny ARAH y P BY BUD I Fr PAS x 的 影响 ; 
ds 代表 栅 格 x 和 y 之 间 的 距离 ;di 代表 威胁 因子 7 
的 影响 范围 ;D, 表 示 生 境 在 胁迫 压力 下 表现 出 退化 
的 程度 ,为 生境 退化 度 指数 ;wr,\B,、5; 依 次 代表 不 
同 威胁 因子 的 权重 、 威 胁 因 子 强 度 、 生 境 抗 干扰 水 
平 . 不 同 生境 对 不 同 威胁 因子 的 相对 敏感 程度 ; 生 
境 质 量 Q, 是 土地 利用 / 敌 盖 类 型 i 中 栅 格 x 的 生境 质 
量 值 ;万 代表 土地 利用 / 黎 盖 类 型 的 生境 适应 比重 ;z 
为 尺度 常数 , 取 值 2.5; 是 半 饱 和 常数 ,在 模型 中 根 


据 栅 格 数据 空间 分 辨 率 (300 m) 来 自 定义 半 饱 和 常 
Beko 计算 出 的 ww 在 0~1 之 间 , 越 接近 1 生境 质量 
越 好 。 

综合 考虑 研究 区 状况 .已 有 学 者 的 相关 人 研究 以 
及 专家 咨询 ,选择 耕地 、 人 造 地 表 、 未 利用 地 作为 胁 
迫 因子 。 根 据 每 个 胁迫 因子 设置 最 大 胁迫 距离 AN 
重 值 和 空间 衰退 类 型 ;用 于 计算 的 7 个 土地 利用 / 窗 
盖 类 型 对 应 胁迫 因子 需 设置 敏感 度 参 数 , 包 括 各 土 
地 利用 / 歼 盖 类 型 的 生境 适宜 度 和 土地 利用 / 歼 盖 类 
型 对 胁迫 因子 的 敏感 度 ”。 参 考 InVEST 模 型 用 户 


表 1 各 景观 类 型 对 生境 威胁 因子 的 敏感 度 
Tab. 1 Sensitivity of landscape types to habitat threat 


factors 
ae 胁迫 因子 
土地 利用 /有 覆盖 类 型 生境 适宜 度 $k JGR GRAIN 

耕地 0.3 0 0.8 0.4 
林地 1 0.6 0.4 0.2 
草地 1 0.8 0.6 0.6 
湿地 0.9 0.6 0.7 0.4 
水 体 0.7 0.5 0.4 0.2 
人 造 地 表 0 0 0 0.1 
未 利用 地 0.6 0.6 0.4 0 


表 2 威胁 因子 及 其 参数 设置 


Tab.2 Threat factors and its parameter setting 


威胁 因子 最 大 胁迫 距离 权重 空间 衰退 类 型 
耕地 6 0.6 线性 
人 造 地 表 9 0.9 Ha 
未 利用 地 4 0.4 线性 


手册 .学 者 们 的 人 研究 成 果 以 及 专家 咨询 结果 进行 参 
数 设置 ,具体 如 表 1 和 表 2 所 示 。 

1.3.2 Sent+MK 44-477 Theil-Sen Median 趋势 
分 析 是 一 种 稳定 且 受 数据 异常 值 影响 较 小 的 计算 
方法 ,为 了 判断 不 同时 间 点 上 收集 到 的 数据 趋势 ， 
通常 结合 具有 正 态 分 布 特征 变量 的 非 参数 Mann- 
Kendall 检验 方法 使 用 ,已 成 为 时 间 序 列 变化 分 析 当 
中 重要 的 人 研究 方法 。 其 对 数据 误差 的 干扰 具有 和 较 
强 的 规避 力 ,不 需要 数据 服从 一 定 的 分 布 , 并 且 为 
显著 性 水 平 检验 提供 了 坚实 的 统计 理论 基础 ,使 结 
果 更 加 科学 可 靠 %71。Theil-Sen Median 趋势 计算 公 
式 为 : 


X,-X, 
Sy= Median[ z | (4) 


j-i 
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式 中 : Sy 表示 生境 质量 指数 变化 趋势 , 当 Sv >0 
时 ,反映 生境 质量 (X) 旦 现 增长 的 趋势 ,反之 则 反映 
生境 质量 (X) 呈 现下 降 的 趋势 。 

Mann-Kendall 的 检验 方法 是 用 来 判断 变化 趋势 
的 显著 性 ,可 适用 于 不 具有 正 态 分 布 特征 变量 的 趋 
势 分 析 。 计 算 公式 如 下 : 


WE X, }, i=1992, 1993, ---, 2020 

S-11 s>0 

/s(S) n-l n 
Z=; 0, S=0 其 中 ， S= >， sgn(X,—X,) 5) 

S+1 <0 7=1 +1 

|s(S) 

eee D2n+5 

sgn(X,—X)= 0, X-X, ,9(8)= n+ 6) 


-1, X,-X, 
AP: X, 和 XX 分 别 表示 栅 格 第 i 年 和 第 j 年 的 生境 
质量 (X) 值 ;n 表示 时 间 序 列 的 长 度 (29 a) ;sgn 为 计 
算 符 号 函数 ;统计 量 Z 的 取 值 范围 在 正 负 无 穷 之 
间 。 在 给 定 显 著 性 水 平 a( 一 般 取 0.05) 下 , 当 
|Z1>uw_ww 时 ,表示 研究 序列 在 a(0.05) 水 平 上 存在 
显著 变化 ”。 

结合 Sen 变化 趋势 程度 和 MK 检验 结果 ,对 生境 
质量 变化 趋势 类 型 分 类 , 当 Sen>0.0005 日 MK>1.96 
时 ,表征 为 生境 质量 明显 上 升 ; 当 Sens-0.0005 H. 
MK<-1.96 时 ,表征 为 生境 质量 明显 降低 ; 当 -0.0005< 
Sen<0.0005 日 -1.96<MK<1.96 时 ,表征 为 生境 质量 稳 
定 不 变 ; 当 Sen>0.0005 有 昌 -1.96<MK<1.96 时 ,表征 为 
生境 质量 轻微 上 升 ; 当 Sen<-0.0005 日 -1.96<MK< 
1.96 时 ,表征 为 生境 质量 轻微 降低 。 
1.3.3 Hurst 指数 发 展 趋势 预测 ” Hurst 指数 用 于 判 
断 时 序数 据 是 否 具有 可 持续 性 或 反 持续 性 ,计算 方 
法 , 详 见 文献 [38]。Hurst>0.5 表 示 生 境 质 量 发 展 趋 
势 与 过 去 趋势 是 一 致 的 , 且 Hurst 值 越 接 近 于 1, 可 
持续 性 越 强 。Hurst<0.5 表 示 生 境 质量 发 展 趋势 与 
过 去 趋势 是 不 一 致 的 , 且 Hurst 值 越 接近 0, 反 持续 
性 越 强 。Hurst=0.5 表示 生境 质量 变化 为 随机 序列 ， 
不 存在 长 期 相关 性 。 

为 了 清晰 揭示 黄河 中 游 生境 质量 的 可 持续 性 
趋势 ,本文 综 合 利 用 Sen 趋 势 .MK 检验 结果 和 Hurst 

站 数 对 生境 质量 发 展 趋 势 进 行 分类 。 当 Sen> 

0.0005 、MK>z1.96 H. Hurst>0.5( Hurst<0.5 ) 时 ,表示 未 
来 生境 质量 发 展 表 现 为 可 持续 性 ( 反 持 续 性 ) 上 升 
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趋势 ; “4 Sens- 0.0005, MK<- 1.96 H. Hurst>0.5 
(Hurst<0.5 ) 时 ,表示 未 来 生境 质量 发 展 表现 为 可 持 
续 性 ( 反 持 续 性 ) 下降 趋 势 ; 当 -0.0005<Sen< 
0.0005 ,-1.96<MK<1.96 H. Hurst>0.5 (Hurst<0.5 ) 时 , 
表示 无 法 确定 未 来 生境 质量 发 展 趋 势 。 


2 结果 与 分 析 


2.1 黄河 中 游 生境 质量 年 际 变化 特征 

从 整个 黄河 中 游 地 区 来 看 ,1992 一 2020 年 生境 
质量 在 0.646~0.657 之 间 发 生 微 弱 的 波动 减少 ,每 
10 a 减少 0.004(P<0.001)。 在 1992 一 2000 年 生境 
质量 在 0.654~0.657 之 间 波 动 减 小 ,每 10 a 减少 
0.003(P<0.001), 在 2001 一 2010 年 生境 质量 减速 程 
度 略 有 加 强 , 每 10 a 减少 0.005 (P<0.001) ;而 在 
2011—2020 年 生境 质量 变化 差异 性 不 显著 (P< 
0.5) ,在 此 10a 间 生境 质量 下 降 趋势 缓解 ,但 波动 变 
化 较 大 (图 2a)。 从 黄河 中 游 退耕 还 林 ( 草 ) 区 分 析 
结果 来 看 (图 2b) ,1992 一 2020 年 生境 质量 呈现 增加 
的 趋势 ,每 10 a 增加 0.174(P<0.001)。 在 1992 一 
2000 年 黄河 中 游 退耕 还 林 ( 草 ) 区 生境 质量 在 
0.304~0.307 之 间 波 动 减 小 ,每 10 a 减少 为 0.010 
(P<0.001) ,而 从 2001 年 该 地 区 退耕 还 林 ( 草 ) 工 程 
的 大 规模 实施 ,2001 一 2010 年 生境 质量 明显 提升 ， 
在 0.308~0.402 之 间 迅 速 增加 ,其 在 每 10 a 增加 
0.122 (P < 0.001) ; 2011—2020 年 生境 质量 显著 提 
升 ,由 0.451 又 增 到 0.966, 增 速达 到 0.656 (P< 
0.05)。 
2.2 黄河 中 游 生境 质量 时 空 变化 趋势 

根据 黄河 中 游 及 退耕 还 林 ( 草 ) 区 的 生境 质量 
年 际 变化 规律 (图 2) ,选取 研究 区 生境 质量 变化 较 
为 明显 的 时 间 段 (2001 一 2020 年 ) ,开展 黄河 中 游 生 
境 质量 时 空 变化 分 析 。 从 2001 年 .2011 年 和 2020 
年 3 个 时 期 来 探究 黄河 中 游 不 同时 期 生境 质量 等 级 
空间 分 布 ,将 生境 质量 指数 按照 0~0.2、.0.2~0.4、0.4~ 
0.6 .0.6~0.8 .0.8~1 划分 为 差 . 较 差 一 般 、 良 好 和 优 5 
TERRES 

黄河 中 游 生境 质量 的 空间 格局 以 优 和 较 差 等 
级 为 主 , 优 等 级 主要 分 布 在 太行 山 Rak 
部 等 区 域 , 较 差 等 级 主要 分 布 在 关中 盆地 ` 渭 河 平 
原 等 区 域 。 然 而 , 随 着 社会 经 济 快速 发 展 , 生 境 质 
量 等 级 为 差 的 范围 在 逐渐 辐射 扩散 (图 3a~ 图 3c)。 


202303.10341v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


460 干旱 区 人 研究 


40 卷 


(a) 


~~ - slope= —0.003, P<0.001 
> -- - slope= —0.005, P<0.001 
- - - slope= —0.001, P<0.5 


年 份 
图 2 1992 一 2020 年 黄河 中 游 (a) 及 


Fig. 2 Change trend of habitat quality in the middle reaches of the Yellow River (a) and the Grain for Green Project 
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~~~ slope=-0.010, P<0.001 
-- - slope=0.122, P<0.001 
~~~ slope=0.656, P<0.05 


年 份 


退耕 还 林 ( 草 ) 区 (b) 生 境 质量 变化 趋势 
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Fig. 3. Spatiotemporal distribution of habitat quality grades in the middle reaches of the Yellow River at different periods 


从 研究 区 生境 质量 等 级 在 不 同时 间 段 的 面积 占 比 
来 看 ,生境 质量 等 级 为 优 的 面积 占 比 均 超过 51%， 
良好 和 一 般 等 级 保持 稳定 ,而 生境 质量 的 较 差 和 差 
等 级 变化 明显 ,其 中 较 差 等 级 面积 占 比 由 47.4% 下 


降 到 45.5% , 面积 减少 约 2% , 差 等 级 面积 则 增加 了 


Ut 2%( Al 3d). 


从 2001—2020 年 ,2001 一 2010 年 ,2011—2020 
年 3 个 时 间 段 ,分析 黄河 中 游 生境 质量 时 空 变 化 趋 


202303.10341v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


3 期 陈 ， 实 等 :退耕 还 林 ( 草 ) 工 程 实施 前 后 黄河 中 游 生 境 质量 时 空 变 化 分 析 461 


势 。 在 整个 黄河 中 游 地 区 ,2001 一 2020 年 有 37203  ” 表 3);2001 一 2010 年 存在 43865 个 像 元 的 生境 质量 
个 像 元 的 生境 质量 发 生 了 变化 ,降低 和 上 升 的 面积 。 发 生 了 变化 ,降低 和 上 升 的 面积 分 别 占 变化 面积 的 
分 别 占 变化 面积 的 81.14% 和 18.86% ,其 中 明显 降低 86.76% 和 13.24% ,其 中 明显 降低 和 明显 上 升 的 面积 
和 明显 上 升 的 面积 分 别 占 80.71% 和 12.19%( 图 4a， 分别 占 84.56% 和 8.48% (图 4b, 表 3) ;2011 一 2020 
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图 4 生境 质量 变化 类 型 空间 分 布 
Fig.4 Spatial distribution of variation types of habitat quality 
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RI 各 时 间 段 黄河 中 游 生境 质量 变化 趋势 类 型 像 元 个 数 及 百分比 
Tab.3 The number and percentage of pixels by type of habitat quality change trend in the middle reaches of the Yellow 


River during each period 


变化 类 型 2001 一 2020 年 2001 一 2010 年 2011 一 2020 年 
变化 像 元 数 /个 百分比 % 变化 像 元 数 /个 百分比 % 变化 像 元 数 /个 百分比 % 
明显 降低 30026 80.71 37097 84.56 25391 57.49 
轻微 降低 161 0.43 964 2.20 931 2:11 
轻微 上 升 2482 6.67 2086 4.76 10588 23.97 
明显 上 升 4534 12.19 3718 8.48 7258 16.43 
总 计 37203 100.00 43865 100.00 44168 100.00 


年 ,研究 区 有 44168 个 像 元 发 生 了 变化 ,生境 质量 呈 
现 上 升 趋势 的 面积 较 2001 一 2010 年 有 所 增加 ,而 呈 
现 降低 趋势 的 面积 较 2001 一 2010 年 有 所 减少 ,分 别 
占 变 化 面积 的 59.60% 和 40.40%, 其 中 轻微 上 升 和 明 
显 上 升 的 面积 大 幅度 增加 ,分 别 占 23.97% 和 16.43%， 
明显 降低 的 面积 又 减 至 57.49%( 图 4c, 表 3)。 

从 空间 格局 来 看 ,2001 一 2010 年 生境 质量 明显 
降低 的 区 域 主要 分 布 在 汾 渭 平 原 快 速 城镇 化 发 展 


3248 个 像 元 的 生境 质量 发 生 了 变化 ,表现 为 生境 质 
量 上 升 ,其 中 轻微 上 升 和 明显 上 升 的 面积 分 别 占 
60.65% 和 39.35%( 图 44, 表 4);2001 一 2010 年 存在 
1452 个 像 元 的 生境 质量 发 生 了 变化 ,主要 表现 为 生 
境 质 量 上 升 ,其 中 轻微 上 升 和 明显 上 升 的 面积 分 别 
占 65.01% 和 34.71%( 图 4e, 表 4);2011 一 2020 年 ,人 研 
究 区 有 4804 个 像 元 发 生 了 变化 ,生境 质量 呈现 上 升 
趋势 的 面积 较 2001 一 2010 年 有 明显 的 增加 ,生境 质 


区 ,如 陕西 西安 市 .陕西 太原 市 片区 ,生境 质量 明显 
上 升 的 区 域 主要 分 布 在 陕西 省 榆林 市 东北 部 、 延 安 
市 东部 局 部 .甘肃 省 天 水 市 西部 局 部 ,以 及 山西 省 
忻州 市 和 太原 市 西部 局 部 等 地 区 ;到 2011—2020 
年 ,生境 质量 明显 降低 的 区 域 有 所 缩减 , 仍 分 布 在 
汾 渭 平 原 快 速 城镇 化 发 展区 ,而 生境 质量 明显 上 升 
的 区 域 呈现 多 点 开花 态势 ,如 陕西 省 延安 市 局 部 、 


量 轻微 上 升 (80.50% ) 比 生境 质量 明显 上 升 的 面积 
占 比 要 大 ,同时 还 存在 很 小 比例 的 明显 降低 像 元 ， 
面积 占 比 仅 为 0.52%( 图 4f, 表 4)。 从 空间 格局 来 
看 ,2001 一 2010 年 生境 质量 上 升 的 区 域 主要 分 布 在 
陕西 省 榆林 市 东北 部 ;2011 一 2020 年 ,生境 质量 上 
升 的 区 域 呈 分 散 式 增加 ,主要 分 布 在 山西 省 忻州 市 
和 晋 中 市 河南 省 济源 市 、 甘 肃 省 庆 阳 市 和 年 西 


甘肃 省 庆 阳 市 和 定西 市 局 部 .山西 省 长 治 市 东部 局 
部 以 及 河南 省 济源 市 局 部 等 地 区 。 总 体 上 ,2001 一 
2020 年 黄河 中 游 地 区 生境 质量 降低 区 域 主要 分 布 
在 快速 城镇 化 发 展 地 区 ,生境 质量 上 升 主 要 分 布 在 
生态 工程 实施 的 重点 地 区 。 

在 黄河 中 游 退 耕 还 林 ( 草 ) 区 ,2001 一 2020 年 有 


市 。 总 体 来 看 ,2001 一 2020 年 退耕 还 林 ( 草 ) 区 生境 
质量 上 升 区 域 主要 分 布 在 陕西 省 榆林 市 东北 部 和 
山西 省 忻州 市 中 西部 。 
2.3 黄河 中 游 生境 质量 变化 可 持续 性 分 析 

黄河 中 游 地 区 生境 质量 的 Hurst 指数 范围 为 
0.318~0.760( 均 值 为 0.685) , 反 持 续 性 (Hurst<0.5) 占 


表 4 各 时 间 段 黄河 中 游 退耕 还 林 ( 草 ) 区 生境 质量 变化 趋势 类 型 像 元 个 数 及 百分比 
Tab.4 The number and percentage of pixels by type of habitat quality change trend in the area of the Grain for Green 


Project in the middle reaches of the Yellow River during each period 


as pate 2001—20204F 2001—20104F 2011—2020 4 

变化 像 元 数 /个 百分比 % 变化 像 元 数 /个 百分比 % 变化 像 元 数 /个 百分比 % 
明显 降低 0 0 4 0.28 25 0.52 
轻微 降低 0 0 0 0 0 0 
轻微 上 升 1970 60.65 944 65.01 3867 80.50 
明显 上 升 1278 39.35 504 34.71 912 18.98 
总 计 3248 100.00 1452 100.00 4804 100.00 
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HEAL 0.88% ,可 持续 性 (Hurst>0.5) 占 比 高 达 99.12%。 
结合 Sen 趋势 和 MK 检验 的 地 理 空间 分 析 来 看 ,黄河 
中 游 生境 质量 在 维持 2001] 一 2020 年 的 变化 趋势 下 ， 
以 持续 性 降低 和 持续 性 上 升 发 展 为 主 ,生境 质量 可 持 
续 性 降低 趋势 区 所 覆盖 区 域 较 大 ,主要 分 布 于 汾 渭 
平原 以 及 城市 集中 发 展区 域 ,生境 质量 可 持续 性 上 
升 趋势 区 主要 位 于 陕西 省 的 榆林 市 北部 .山西 省 忻 
州 市 西部 局 部 ,甘肃 省 的 天 水 市 和 平凉 市 局 部 (图 
Sa)。 
在 黄河 中 游 退耕 还 林 ( 草 ) 区 ,生境 质量 的 Hunrst 
旨 数 范围 为 0.483~0.760( 均 值 为 0.644) ,可 持续 性 
(Hurst>0.5) 占 比 达到 了 99.09%。 未 来 年 份 黄河 中 
游 退耕 还 林 ( 草 ) 区 生境 质量 将 维持 2001 一 2020 年 
的 变化 趋势 ,体现 在 生境 质量 可 持续 性 上 升 ,主要 
分 布 在 陕西 省 (榆林 市 )、 山 西 省 (忻州 市 ) 以 及 内 蒙 
古 自治 区 ( 乌 兰 察 布 市 .呼和浩特 市 和 鄂尔多斯 市 ) 
三 地 交界 片区 (图 5b)。 


3 讨论 


在 InVEST 模 型 估算 中 ,草地 和 林地 的 生境 适宜 
度 高 于 人 造 地 表 和 耕地 , 随 着 退耕 还 林 ( 草 ) 生 态 工 
程 实施 ,耕地 向 林地 或 草地 转换 的 面积 逐年 增加 ， 
黄河 中 游 生境 质量 年 际 变化 伴随 着 耕地 向 林地 和 
草地 面积 变化 而 变化 。 从 黄河 中 游 区 2001 一 2010 
年 和 2011 一 2020 年 两 个 时 间 有 段 的 土地 利用 / 窗 羔 类 
型 转移 情况 来 看 ,2001 一 2010 年 耕地 仅 有 3836 个 像 
元 向 林地 (621 个 像 元 ) 和 草地 (3215 个 像 元 ) 转 移 
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(图 6a) ,以 致 退耕 还 林 ( 草 ) 区 2001 一 2010 年 生境 质 
量 上 升 缓慢 ,而 在 2011 一 2020 年 耕地 向 林地 (3175 
个 像 元 ) 和 草地 (18430 个 像 元 ) 转 移 的 像 元 增加 到 
21605 个 (图 6b) ,在 此 阶段 黄河 流域 退耕 还 林 ( 草 ) 
区 的 生境 质量 显著 上 升 。 

考虑 到 黄河 中 游 生境 质量 年 际 变化 规律 ,以 及 
退耕 还 林 ( 草 ) 生 态 工 程 全面 实 施 的 时 间 节 点 ”, 划 
分 了 2001—2020 年 .2001 一 2010 年 和 2011 一 2020 
年 时 间断 面 ,从 整体 (黄河 中 游 地 区 ) 和 局 部 (黄河 
中 游 退耕 还 林 ( 草 ) 区 ) 两 个 尺度 ,在 退耕 还 林 ( 草 ) 
工程 实施 前 后 ,探究 生境 质量 时 空 变 化 特征 。 在 整 
体 上 ,黄河 中 游 地 区 生境 质量 总 体 较 好 ,但 在 年 际 
变化 上 表现 为 轻微 的 波动 下 降 , 这 与 杨 洁 等 "研究 
黄河 流域 的 平均 生境 质量 指数 0.631 ,整体 呈现 小 幅 
下 降 趋 势 的 结果 相近 。 在 局 部 上 ,黄河 中 游 退耕 还 
林 ( 草 ) 区 生境 质量 在 年 际 变化 上 呈现 明显 上 升 的 
态势 ,这 与 Guo 等 探究 呼伦贝尔 一 包头 一 鄂 尔 多 
斯 一 榆林 城市 群 和 太原 城市 群 的 生态 系统 服务 价 
值 所 呈现 的 上 升 趋势 较为 一 致 。 

从 不 同 空间 尺度 研究 均 可 发 现 , 由 耕地 转 为 林 
地 、 草 地 的 区 域 ,其 生境 质量 上 升 发 展 的 可 持续 性 
较 好 ,黄河 中 游 区 计划 和 正在 实施 的 退耕 还 林 ( 草 ) 
生态 工程 地 区 具备 生境 质量 上 升 潜力 ” ,生境 质量 
上 升 地 区 主要 位 于 榆林 、 朔 州 等 生态 工程 实施 的 重 
不 地 区 ,这 些 地 区 本 身 属于 经 济 欠 发 达 的 山区 , 受 
人 类 活动 干扰 较 小 ,大 面积 的 押 荒 土地 增加 了 植被 
ARU ,退耕 还 林 ( 草 ) 生态 文明 建设 等 政策 措 
施 提 高 了 生境 质量 汪 。 在 研究 区 城镇 原 建成 区 的 
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图 5 2001 一 2020 年 生境 质量 时 空 变化 趋势 
Fig.5 Spatiotemporal change trend of habitat quality from 2001 to 2020 


202303.10341v1 


chinaXiv 


464 FẸ 


ChinaXiv 合 作 期 刊 


(a) 2001 一 2010 年 


(b) 2011 一 2020 年 
wal 
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Fig. 6 Land use and cover type transfer in the middle reaches of the Yellow River in different periods 


生境 质量 表现 为 最 低 等 级 ,但 随 城 镇 化 发 展 , 原 城 
镇 周边 区 域 土地 利用 /覆盖 类 型 不 断 向 城镇 用 地 转 
化 3, 生 境 质 量 下 降 地 区 主要 分 布 在 西安 .太原 等 
城市 及 市 区 周边 , 究 其 原因 是 这 些 地 区 交通 干线 密 
集 、 人 口 密度 大 ,在 经 济 发 展 过 程 中 土地 利用 类 
型 不 断 向 城镇 用 地 转化 ,高 强度 和 不 合理 的 人 类 活 
动 加 剧 ,引发 生境 质量 退化 ,这 些 地 区 未 来 需要 
重点 关注 。 

黄河 中 游 地 区 正 坚 定 不 移 走 生态 优先 .绿色 发 
展 之 路 ,全面 推动 退耕 还 林 ( 草 ) 等 生态 工程 的 实 
施 。 结 合 研 究 结 果 提 出 以 下 几 点 建议 :(1) 生境 质 
量 在 正 向 和 负 向 上 极 易 受 人 类 活动 的 干扰 , 需 针 对 
水 源 涵养 区 .水土 保 持 区 、 湿 地 保护 区 以 及 生态 廊 
道 , 合 理 规划 布局 ,使 经 济 社会 发 展 与 资源 承载 力 
协调 发 展 ;(2) 退耕 还 林 ( 草 ) 能 够 有 效 提升 生境 质 
量 , 治 理 与 保护 并 重 切 实 做 好 封 山 禁 牧 ,实施 宜 林 
则 林 、 宜 草 则 草 的 绿化 工程 ;(3) 城市 及 周边 区 域 生 
境 质 量 降低 ,主要 是 由 于 在 开发 利用 自然 资源 时 ， 
过 于 注重 短期 的 经 济 价值 ,忽视 了 生境 质量 的 保 
护 ,生境 质量 的 提升 与 植被 状况 紧密 相关 ,通过 绿 
化 等 措施 ,增加 林 草 植被 覆盖 度 ,可 有 效 提 高 区 域 
生境 质量 ”。 

本 研究 在 数据 上 使 用 了 时 间 分 辩 率 较 高 的 
CCI-LC 逐年 土地 覆盖 数据 ,但 其 空间 分 辨 紊 有待 提 
升 ,今后 可 在 Google Earth Engine (GEE) 平 台 , 基 于 
Landsat 影像 解 译 的 30 m 土地 覆盖 数据 估算 生境 质 
量 , 提 高 生境 质量 评估 空间 尺度 。 本 研究 采用 的 In- 


VEST 模 型 估算 生境 质量 时 , 仪 考虑 了 威胁 因子 对 
生境 的 累积 胁迫 作用 ,而 威胁 因子 间 的 交互 胁迫 作 
用 却 未 考虑 其 中 ,各 景观 类 型 对 生境 威胁 因子 的 敏 
感度 存在 人 为 主观 干扰 ,人 研究 区 内 部 和 外 部 的 威胁 
因子 没有 得 到 很 好 的 兼顾 。29 a 来 黄河 流域 中 游 地 
区 生境 质量 整体 较为 稳定 ,退耕 还 林 ( 草 ) 工 程 实施 
前 ,生境 质量 的 时 空 变化 强度 较 小 ,但 随 着 生态 工 
程 实施 和 社会 经 济 发 展 , 人 类 活动 持续 增强 ,生境 
质量 在 城镇 及 周边 地 区 表现 出 聚集 扩张 式 的 退 
化 。 未 来 可 将 黄河 中 游 生境 质量 时 空 变 化 影响 机 
制作 为 研究 重点 ,定量 探测 其 驱动 因素 ,以 期 科学 
有 效 的 提出 生态 环境 保护 和 恢复 的 措施 。 


4 结 论 


本 文 基于 黄河 中 游 CCI-LC 逐年 十 地 覆盖 数据 
集 , 依 托 InVEST 模 型 计算 逐年 生境 质量 ,运用 Sent+ 
MK 检验 .Hurst 指 数 分 析 方 法 ,在 退耕 还 林 ( 草 ) 前 后 
探究 了 1992 一 2020 年 黄河 中 游 生境 质量 的 时 空 变 
化 特征 和 发 展 趋势 ,研究 成 果 将 为 黄河 中 游 地 区 开 
展 生态 环境 保护 恢复 提供 科学 指导 。 具 体 结论 
如 下 : 

(1) 近 29 a 来 ,黄河 中 游 生 境 质 量 呈 现 微弱 波 
动 减 小 ,而 在 其 退耕 还 林 ( 草 ) 区 ,生境 质量 呈现 明 
显 增加 的 趋势 ,其 中 在 退耕 还 林 ( 草 ) 工 程 大 规模 实 
施 后 ,生境 质量 增加 迅速 。 

(2) 黄河 中 游 生境 质量 时 空 分 布 趋势 变化 的 降 
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低 类 型 所 占 面 积 比 上 升 类 型 要 大 得 多 ,但 随 着 退耕 
还 林 ( 草 ) 工 程 的 逐步 实施 ,2011 一 2020 年 较 2001 一 
2010 年 生境 质量 呈现 上 升 趋势 的 面积 有 所 增加 。 
尤其 是 在 退耕 还 林 ( 草 ) 区 ,2011 一 2020 年 较 2001 一 
2010 年 生境 质量 上 升 趋势 的 面积 明显 增加 ,轻微 上 
升 和 明显 上 升 类 型 占 主导 。 


(3) 未 来 黄河 中 游 生境 质量 以 持续 性 降低 和 持 


续 性 上 升 发 展 为 主 , 在 退耕 还 林 ( 草 ) 区 生境 质量 以 
可 持续 性 上 升 为 主 , 主 要 分 布 在 榆林 市 .忻州 市 S 


aL 


察 布 市 .呼和浩特 市 和 鄂尔多斯 市 等 地 区 。 


一 
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Spatiotemporal changes of habitat quality before and after the implementation 
of Grain for Green Project in the middle reaches of the Yellow River 


CHEN Shi, HUANG Yinlan, JIN Yunxiang’ 
(1. School of Geography and Planning, Application Research Center of Remote Sensing for Natural Resources, 
Chizhou University, Chizhou 247000, Anhui, China; 2. Institute of Agricultural Resources and Regional Planning, 
Chinese Academy of Agricultural Sciences, Beijing 100081, China) 


Abstract: It is of great practical significance for the ecological protection, rehabilitation, and high- quality 
development of the middle reaches of the Yellow River to analyze the spatial and temporal characteristics of 
habitat quality and its development trend. On the basis of the land cover data of the European Aeronautics and 
Space Administration from 1992 to 2020, the InVEST model was used to estimate the habitat quality. Sen+ MK 
test and Hurst index analysis methods were applied to explore the spatial and temporal changes and development 
trends of the habitat quality in the middle reaches of the Yellow River before and after the implementation of the 
Grain for Green Project. Results showed that the habitat quality in the middle reaches of the Yellow River was 
good overall in recent 29 years, with the average habitat quality index of 0.652. However, a slight downward 
trend was observed (a drop of 0.004 per decade). Meanwhile, the habitat quality in the Grain for Green area 
showed an upward trend (an increase of 0.174 per decade), of which the most significant reached 0.656 recorded 
from 2011 to 2020. The spatial and temporal changes of habitat quality in the middle reaches of the Yellow River 
from 2001 to 2020 mainly showed an obvious decrease trend. In terms of the development of different time 
periods, the area with good habitat quality in 2011-2020 increased compared with that in 2001-2010. This 
finding was reflected in the improvement of habitat quality in the Grain for Green area. In the future, the temporal 
and spatial changes of habitat quality in the middle reaches of the Yellow River will be dominated by obvious 
decrease and obvious increase trends, respectively, and the sustainability will be strong. The significantly reduced 
areas will be mainly distributed in the Fenwei plain and the urban centralized development areas, and the 
significantly increased areas will be mainly located in the border areas of Shanxi Province, Shaanxi Province and 
Inner Mongolia Autonomous Region. 

Keywords: habitat quality; INVEST model; the Grain for Green Project; CCI-LC; the Yellow River Basin 


